Refraction of Gluon jet in QGP 
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We introduce photon-like refraction for gluon jet in quark gluon plasma created in high energy 
nucleus-nucleus collisions. A quark jet goes straightly while a gluon jet is assumed to bend in 
the non-homogeneous medium. Combining refraction and energy loss, we calculate the away-side 
dihadron A(/> distribution in central Au + Au collisions at RHIC energy. The dihadrons from quark- 
quark pair center around A(j> — it, while the dihadrons from quark-gluon and gluon-gluon pair are 
deflected from A<j!> = n. The resulting distribution shows three peaks at intermediate pr since both 
quark-quark and quark-gluon are important. However we only observe one center peak at high 
Pt, where the quark-quark contribution dominates. These observation qualitatively agrees with the 
experiment data. 

PACS numbers: 12.38.Mh, 24.85.+p; 25.75.-q 



A strongly couple partonic matter (sQGP) is now be- 
lieved to exist in high energy nucleus-nucleus collisions 
with an extended volume and time. The high trans- 
verse momentum px partonic jets traversing the hot 
dense matter lose a significant amount of their energy 
via induced gluon radiations, a phenomena known as jet- 
quenching PJ. Such energy loss is predicted to lead to 
strong suppression of both single- and correlated away- 
side dihadron spectra at highp r @, H, I3| i consistent with 
experimental findings [ESQ- 

The suppression factor of the leading hadrons from 
jet fragmentation depends on the total parton energy 
loss which in turn is related to the jet propagation path 
weighted with the gluon density p g 0, 0] , 
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for a parton produced at a transverse position r and trav- 
eling along the direction n. The energy dependence of the 
energy loss is parameterized as 
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from the numerical results in Ref. Q in which thermal 
gluon absorption is also taken into account. The param- 
eter eo should be proportional to the initial gluon den- 
sity. It was shown that the single hadrons at high pt 
are dominated by the jets emitted from surface layer of 
the overlap. However, the away-side dihadrons contains 
contribution from both tangentially emitted or punch- 
through dijets. A simultaneous fit to the single and di- 
hadron data constrains the energy loss parameter within 
a narrow range: eo = 1.6 — 2.1 GeV/fm 

Recently the azimuthal angle A(f> correlations are pre- 
sented for charged hadrons from dijets for 0.4 < pt < 
15 GeV/c in Au+Au collisions at ^/snn = 200 GeV 
by PHENIX group. Especially, with increasing px, the 
away-side distribution evolves from a broad to a concave 



shape, then to a convex shape Q. Comparisons to p + p 
data suggest that the away-side can be divided into a par- 
tially suppressed "head" region centered at A0 ~ n, and 
an enhanced "shoulder" region centered at A<fi ~ 7r± 1.1. 
These away-side results are then used to constrain various 
proposed mechanisms, such as mach-cone Q , large angle 
gluon radiation II CI I 111 , Cherenkov gluon radiation [121 ] 
or deflected jets [ljj, Il4 |. The goal of current letter is to 
propose gluon refraction as yet another mechanism for 
the away-side behavior, based on a NLO pQCD parton 
model with jet quenching in A + A collisions. 

For an optical ray propagating in non-homogeneous 
medium with a refractive index distribution n(r). Its 
pathway in the medium should satisfies the relation 
(5 / n(r)ds = 0, which yields the following equation, 
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ds 



n(r) 



Vra(r). 



(3) 



Generally, an optical ray is refracted in a nonho- 
mogeneous medium, but electrons do not. This moti- 
vate us to assume an energetic gluon jet is refracted 
while a quark jet always goes straightly in the non- 
homogeneous medium. On a microscopic level, lattice 
data suggests [l5| that there could be massive bound 
states in strongly-couple QGP around the critical tem- 
perature. These bound states could serve as colored scat- 
tering centers for gluons, and lead to refraction. 

For a photon with energy E, the refractive index n(E) 
is related to the density of the scattering centers in 
medium N, n(E) = 1 + ^-FJE), where F(E) is the 
forward scattering amplitude [161 ]. Similarly, we assume 
the refractive index for a gluon jet is related to the initial 
gluon density of the medium, 



n(E, t, b, r) = 1 + k(E)p s (t, b, r), 



(4) 



where the k(E) factor should be related to the disper- 
sive power of a gluon jet with a given energy, similar to 
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FIG. 1: The trajectories of four types of jet pairs for quark- 
quark (a), quark-gluon (b), gluon-quark (c) and gluon-gluon 
(d), respectively, back-to-back created in central nucleus- 
nucleus collisions. A and B represent the direction of trigger 
and away-side, respectively. 

2nF(E)/ E in Ref.pjj]. Here for a simple calculation, wc 
choose a fixed value of n(E) by fitting data. 

Neglecting transverse expansion, the gluon density dis- 
tribution in a ID expanding medium in A + A collisions 
at impact-parameter b is assumed to be proportional to 
the number of binary collisions, 

p 9 (r,b,r) = — —^Lt A (\r + nT\) 

X t A (|b-(r + nr)|). (5) 

A quark jet produced at r always hold its initial di- 
rection n while a gluon jet produced at r is deflected 
from its initial direction in the medium before fragment- 
ing into hadrons outside the overlap region, tyt(f) = 
2rrB 2 \/ ^ ~ f 2 / R 2 is the nuclear thickness function in a 
hard-sphere geometry model. We use a set of parame- 
ters same as those in Ref. fio = 1-5 GcV, eo^o = 0.5 
GeV and r = 0.2 fm/c in Eq. © © © . The total 
energy loss of a quark or a gluon jet is all given by 
Eq. ©©©, except that the gluon jet follows a curved 
path given by Eq. © © . In the most central Au + Au 
collision, the refractive index can be simply written as 
n(r) = 1 + 2k(1 - r 2 /R 2 ). When we choose k = 0.25 
and energy loss parameter eo = 1.7 GeV/fm, numerical 
results for single and dihadron spectra [l7| can fit data 
well. 

Figure © shows the trajectories for the four possi- 
ble hard-scattered jet pairs, i.e. quark-quark (a), quark- 
gluon (b), gluon-quark (c) and gluon-gluon (d). The 
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FIG. 2: (color online). The single hadron (dashing curve) and 
dihadron (solid curve) cross sections only contributed respec- 
tively from single gluon jets and qg/gq jet pairs as a function 
of 8 of gluon jets deflected from its initial direction when the 
jets arrive at the system surface. 
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FIG. 3: (color online). The dihadron yield per unit area as a 
function of initial spatial place r of surviving and fragmented 
jet pairs in the transverse plane in most central Au + Au 
collisions. 



refraction results in a deflection of the gluon, thus the 
quark-gluon and gluon-gluon jet pairs are no longer back- 
to-back. Consequently, this a acoplanarity propagates to 
the away-side hadron pairs after fragmentation. In ad- 
dition, the gluon deflects towards higher density region, 
leading to a longer path and more energy loss for a gluon 
jet, so its contribution to observed hadron pairs is re- 
duced. Figure © shows the distribution of the deflec- 
tion angle 0, for single hadron from gluons and hadron 
pairs from qg/gq di-jets. The single hadron distribution 
peaks at 8 = 0, but has a long tail, this is consistent 
with the surface emission picture, where the observed jets 
traverse short path in the medium, thus have small de- 
flection angles. However, the dihadron from quark-gluon 
dijet peaks around 6 C = 0.72, implying that gluon jet se- 
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FIG. 4: (color online). Per-trigger yield versus A<f> for various 
trigger and partner pr (j>t ®Pt), arranged by increasing pair 
momentum (sum of and Pt), in p + p and most central 
Au + Au collisions. The Data are from Ref. 0. 



FIG. 5: (color online). Comparisons of R qq /(qg+ 3 g) and Rhs 
as functions of the transverse momentum of the associated 
hadron with increasing trigger bins in — 20% Au + Au col- 
lisions at 200 GeV. The Data are from Ref. 0. 



lected by the hadron pair has to go through on average 
a finite medium, thus suffers a substantial deflection. 

Figure shows the distribution of the origin the jet 
pairs contributing to the observed dihadrons as a func- 
tion of the initial spatial place r in the transverse plane 
for b=0fm. Since the energy loss of a gluon jet is about 2 
times of that of a quark jet, gluon-gluon contributions are 
greatly suppressed and can be ignored in the most cen- 
tral Au + Au collisions. The quark-gluon contributions 
come dominatingly from the outer layer of medium. The 
dihadron yields in Figure @ © have been averaged over 
the azimuthal angle of hadron pairs. Because a gluon 
jet has a biggest deflected angle if created with a moving 
direction perpendicular to the gradient of the index of 
refraction, the descriptions for q — g contributions in Fig- 
ure ©© are consistent to the tangential surface emis- 
sion bias. 

The NLO framework of Ref. [TtJ is used to calculate 
the away-side dihadron multiplicity separately for q — q, 
q — g and g — g process. However, the azimuthal distribu- 
tion of dihadron in a NLO pQCD is divergent at A</> = it 
since there are large logarithmic high order corrections 
(Sudakov logs) which need be resummed [H|]. Such diver- 
gence is generally treated by smearing A<j> to match the 
p+p shape but keep the integral over A</> fixed. In this 
paper, we follow the same approach for the case of p + p 
collisions, and for quark-quark contributions in Au + Au 
collisions by simply smearing around A<p = n. Because 
of gluon refraction in Au + Au collisions, the smearing 
technique isn't needed for quark-gluon and gluon-gluon 
contributions, as shown in Figure © where the Gaussian 
curve is directly calculated out without any smearing. 

The final per-trigger yield A(f> distribution are shown 
in Figure ^ . The calculation is compared with PHENIX 
data in several trigger "A" and partner "B" pt {Pt®Pt) 
bins. Our calculation clearly shows one center peak 



around A<j) — tt and two side peaks around A<p = 7r±0.72 
at the away-side. The center peak ( "head" region) is do- 
mianted by the quark-quark jet pairs, while the two side 
peaks ( "shoulder" region) are dominated by quark-gluon 
and gluon-gluon dijets. However for the 5 — 10 <E> 5 — 10 
GeV/c bin, the side peak amplitude decreases relative 
to the center peak. This is because gluon jet contribu- 
tions to hadrons decrease with increasing hadron px, thus 
the dihadron yield is dominated by the quark-quark jet 
pairs. This effect can also be seen in Figure (|3J|. Our 
side-peak location of 0.72 is still smaller than the experi- 
mental value of 1.1, suggesting there might be additional 
deflection mechanism. 

We can also quantify the relative importance of the 
side peaks and center peak using the Rhs introduced by 
PHENIX [3]. In our case, it simply reflects the relative 
contribution from quark-gluon process and quark-quark 
process (neglecting contributions from gluon-gluon), i.e, 

/ dA<j){dN AB /dA4>)ii 
Rhs ~ JdA0(dN AB /dAq>)ig+M ~ Rqq/qa (6) 

Rqq/q g is insensitive to k(E) and n(E, r, b, r) accord- 
ing to Eq. Q. Figure ([5]) shows the comparison of the 
Rqq/qg with the experimental data. The excellent agree- 
ment suggests that gluon deflection is reasonable. And 
the Rhs largely reflects the ratio of the surviving quark- 
gluon to quark-quark in central collisions. 

We have repeated similar calculations for LHC energy. 
Since both gluon-gluon and quark-gluon dominates over 
quark-quark up to pt ~ 100 GeV/c, we expect Rhs < 
1. Since gluon-gluon has stronger deflection than quark- 
gluon, there might even exist multiple side peaks [17] . 

In summary, we argue that the refraction for gluon 
jet might be a possible mechanism for the rich and pt 
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dependent away-side shape in central Au+Au collisions 
observed at RHIC. The gluon refraction is implemented 
in the NLO pQCD energy loss framework to obtain the 
away-side A0 distribution. The away-side shows a center 
peak in the head region and two side-peaks in the shoul- 
der region. The former is dominated by the quark-quark 
dijets while the latter is dominated by the qg/gq dijets. 
The different proportions of contributions from qq, qg/gq 
and gg pairs determine whether there are one, or three 
peaks in dihadron A</> distribution. Our calculation qual- 
itatively describes the experimental findings, however the 
refraction angle under-predicts the split angle seen in the 
data. This might be due to simplification in our assump- 
tion or additional mechanisms need to be included. If 
gluon jets are indeed refracted, one might expected dif- 



ferent shape at LHC at the same p?. Also since the high 
Pt (anti-)proton production is dominated by hard gluon 
fragmentation, the azimuthal angle distribution of the 
high (anti-)proton-(anti-)proton pairs might have dif- 
ferent shape from that of the charged hadron pairs mea- 
sured at RHIC. These studies will be done in our future 
works. 
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